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Abstract: When the state of a quantum system belongs to a N -dimensional 
Hilbert space, with N the power of a prime number, it is possible to associate to 
the system a finite field (Galois field) with N elements. In this paper, we introduce 
generalized Bell states that can be intrinsically expressed in terms of the field oper- 
ations. These Bell states are in one to one correspondence with the elements of 
the generalised Pauli group or Heisenberg- Weyl group. This group consists of dis- 
crete displacement operators and provides a discrete realisation of the Weyl function. 
Thanks to the properties of generalised Bell states and of quadratic extensions of 
finite fields, we derive a particular solution for the Mean King's problem. This solu- 
tion is in turn shown to be in one to one correspondence with a set of N'^ self-adjoint 
operators that provides a discrete realisation of the Wigner quasi- distribution. 



1 Introduction 

We showed previously, ^ in the framework of quantum cloning, that in dimension 
= 4 different classes of Bell states can be defined, that are associated to different 
groups of permutations of the basis states. These Bell states were shown to be in 
one to one correspondance with a commutative group (generalised addition) of N 
elements. More recently, we showed thatin dimension 4 there exists a multiplica- 
tion that together with the addition forms a finite field or Galois field and that the 
generalised Bell states can be defined intrinsically in terms of the field properties. 

We also showed that whenever the dimension of the Hilbert space is a prime 
power (A^ = with p a prime and m a positive integer), these properties can be 
generalised. This is due to the fact that finite fields (sometimes called division rings) 
with elements only exist when A^ is a prime power. 

Besides, it is knownP', that when the dimension is a prime power there exists 
a maximal set of A^ + 1 mutually unbiased El bases or MUB's (two orthonormal 
bases of a A^ dimensional Hilbert space are said to be mutually unbiased if whenever 
we choose one state in the first basis, and a second state in the second basis, the 
modulus squared of their in-product is equal to 1/A^). We showed in Ref. how to 
(re) derive in a straightforward manner an expression for the states of these bases, in 
agreement with the expressions of Wootters and Fields ^ and with the approach of 
Bandyopadhiay 4j et al. . This expression can be intrinsically formulated in terms 
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of the field operations (addition and multiplication between elements) and of the 
additive characters of the field (pth root of unity when = p"^). 

The existence of + 1 MUB's is directly related to the so-called Mean King's 
problem ^ ^Oj which can be formulated as follows, in the qubit case (A^ = 2): 

A Mean King challenges a physicist, Alice, who got stranded on the 
remote island ruled by the King, to prepare a spinl/2 atom in any state 
of her choosing and to perform a control measurement of her liKing. 
Between her preparation and her measurement, the King's men deter- 
mine the value of either ax, or oz- Only after she completed the 
control measurement, the physicist is told which spin component has 
been measured, and she must then state the result of that intermediate 
measurement correctly. How does she do it? 

In other words, how is it possible to ascertain the spin values of the spin observ- 
ables along the X, Y and Z directions? The corresponding eigenbases are MUB's, 
so that the problems consists of ascertaining the values of 2 + 1 non degenerate 
observables of which the eigenbases are mutually unbiased relatively to each other. 

In prime power dimension (A^ = p™), the Mean King's problems consists of 
finding a way to ascertain the values of A^ + 1 non-commuting and non-degenerate 
observables that are diagonal in the A^ + 1 mutually unbiased bases. 

We shall show in the present paper a way to generalise to prime power dimensions 
the solutions that can be found in the litterature for the Mean King's problem in 
dimension 2 and in prime dimensions, making use of the properties of the generalised 
Pauli group, of the generalised Bell states, and of their transformation law. This 
solution is shown to be a special case of the general solution obtained by AravindjlUj. 

An appealing property of our particular solution is that it provides a discrete 
counterpart of Wigner's distribution that is valid in even and odd prime power 
dimensions as well. We shall also show in appendix how it is possible to generalise 
certain of these results in arbitrary odd dimensions, provided we reformulate the 
Mean King's problem in a slightly different manner. 



2 Preliminary concepts 

Whenever the dimension A of a Hilbert space is a prime power (A = p*"), with p 
a prime number, and m a positive integer, it is possible to associate to the Hilbert 
space a finite field with A elements^. 

field is a set with a multiplication and an addition operation which satisfy the usual rules, 
associativity and commutativity of both operations, the distributive law, existence of an additive 
identity and a multiplicative identity 1, additive inverses, and multiplicative inverses for every 
element, excepted. 
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We shall in the following label the elements of the field by a muple of integer 
numbers {iojii, ■■■,im-i) running from to p — 1. This muple {io,ii, ...,im-i) is in 
turn in one to one correspondence with an integer number i, < i < N — 1) that 
we define by its p-ary expansion as follows 

In what follows we shall make no difference between elements of the fields and 
their integer counterpart, to the contrary of usual conventions. 

Any finite field is characterized by two operations, a multiplication and an ad- 
dition, that we shall denote Qq and ®g respectively. 

The index G refers to Evariste Galois and is introduced in order not to confuse 
these operations with the usual (complex) multiplication and addition for which no 

index is written. 

It is not a simple problem to find the multiplication table of finite fields, but 
mathematicians solved the problems and such tables are available "on line" . 

Besides, we can assign an integer label to the elements of the field in such 
a way that when the elements of the field are in one-to-one correspondence with 
muples (as explained before), the addition is equivalent with the addition modulo 
p componentwise. Actually this procedure corresponds to what is called a choice 
of basis for the field. It imposes certain contraints on the identification between 
integers and elements of the field. For instance, we must assume in order to avoid 
inconsistencies that corresponds to the neutral element for addition. It is a natural 
choice to associate 1 to the neutral element for multiplication; then , the m — 1 
remaining powers of p are assigned to elements of the field that are independent in 
the same way that elements of a vectorial space are independent. 

As a consequence of fact that the addition is equivalent to the addition modulo p 
componentwise, one can show that whenever we add an element of the field ptimes 
with itself we obtain 0. In particular, the characteristics of the field, which is 
the smallest number of times that we must add the element 1 (neutral for the 
multiplication) with itself before we obtain (neutral for the addition), is always 
equal to a prime number {p when N = p^). This last property is valid for any finite 
field and ultimately explains why the number of elements of a finite field necessarily 
equals a power of prime. 

Let us denote 7g the pth root of unity: 7g = e^-'^'^l'^ . Exponentiating 7g with 
elements g of the field (with the usual rules for exponientation) , we obtain complex 
phasors of the type 7^ (0 < < A?^ — 1). Such phasors can take p different values. 
They can be considered as a p-valued generalisation of the (binary) parity operation 
^i.{2-K/2).g ^Yi^ox corresponds to the qubit case in the sense that the phasor 7^ (0 < < 
A" — 1) only depends on the value of the first component Qq of the p-ary expansion of 
g. go is nothing else than the remainder of g after division by p, when the division 
by p is taken in the usual sense. 
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Exponentiating 7^, the pth root of unity, by this number we obtain a complex 
phase which is also equal to ■jg exponentiated by the rest of this number after 
division by p. 

It is worth noting that all the finite fields are equivalent, up to a relabelling, and 
that our choice of labelling (and of basis) is dictated by requirements of simplicity 
and convenience. Other approaches are also valid, for instance the ones that imply 
the field theoretical Trace. 3, 4, llj Each approach presents its own advantages 
and disadvantages regarding practical applications. With our conventions, many 
computations are formally similar in prime and prime power dimensions, and certain 
expressions are simpler than in other approaches. The price to pay is a loss of 
generality, in the sense that our conventions can be shown to be equivalent to the 
ones made in the Trace approach provided we perform a special choice of basis in 
the latter approach. !?] 

Making use of the fact that the addition is the addition modulo p, component- 
wise, we can derive the following identity which will appear to be very useful in the 
following: 

lh-lh = 7^ ^^'^ (1) 

Indeed, -fh ■ ih = 1g^'^ = 7g ^'"^ = Ig®"'^' = iT"'^ ■ This relation expresses that 
pth roots of unity are additive characters of the Galois field. [T^ 

The following identity also appears to play a fundamental role in our approach: 

E 7r = A^^.o (2) 



The proof is given in Ref . [7j . 

It is important to note, in order to avoid confusions, that different types of 
operations are present at this level: the internal field operations are labelled by 
the lower index G. They must not be confused with the modulo operations. It 
is worth noting that the property Y^p=Q ^ip®i) = NSg^ is true for the modulo N 
multiplication as well, but 7 must be taken to be equal to the A^th root of unity 
in this case. In prime dimensions 7g is the A^th root of unity and the Galois and 
modulo N operations coincide. In prime power but non-prime dimensions, this is 
no longer true. Because of this, certain applications are also valid even when the 
basic operations, addition and multiplication, do not form a field. This is the case for 
instance with the transformation law of the Bell states that can be derived in all odd 
dimensions, provided we operate with the usual (modulo) operations. In appendix 
1 we show that this property opens the way to a reformulation of the Mean King's 
problem that is valid when the dimension is an odd number, not necessarily equal 
to the power of an odd prime. 
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3 MUB's in prime power dimensions. 



In Ref.fT], we obtain explicitly the following expression for MUB's in prime power 
dimensions: 

1 N-l 

|el) = ^T. 7r^^-^(7j5^-'^^^^^^^¥|e;) (3) 

where le^) represents the fcth basis state of the zth MUB [k runs from to A^ — 1 and 
i from 1 to A^) expressed in terms of the states of the computational basis |e°). The 

N bases so-defined and the computational basis are mutually unbiased relatively to 
each other. 

At this level there remains some ambiguity because a square root is always de- 
fined up to a global minus sign. We showed in Ref. [7j how to rise this ambiguity and 

obtained the following determinations for the square root factor 

• in odd prime power dimensions, (^^(^~^)®<^'?®<^'?))5 = ^(^~i)®Gg0Gg)/G2 ^^]-^gj^g 
represents division in the Galois field). The corresponding expression for the 
MUB's is thus, in odd prime power dimensions, 

■1 N-l 

\el) = 4^T. 7g«^««^7r'^'^"''^"^^/"')|e;), (4) 

• in even prime power dimensions, 

C e(O-l)0G<70G'7)a^/ e((i-l)0Gg0Gg)a ^ n"*"^ ^(i-l)0G2"0G2"^O-l)®G2"0G2"' 

It: 

where i represents e~, and where the coefficients /„ are unambiguously defined 
by the p-ary (here binary) expansion of /: / = Y^=q ^n2", while n' is the 
smallest integer strictly larger than n such that /„/ 7^ 0, if it exists, otherwise; 
here 7g = —1- The corresponding expression for the MUB's is thus, in even 
prime power dimensions (2™), 

1 N-l I 

V A/ ,=0 

where q = Y^^=q Qn'^^, while n' is the smallest integer strictly larger than n 
such that Qn' 7^ 0, if it exists, otherwise. 

Although we proved in Ref. (Tj by direct computation that these bases are mu- 
tually unbiased, in the present paper we shall rather follow the approach of Ref. |3] 
which is more general. 
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In the reference P| it is shown that when there exists a maximal commuting ba- 
sis of orthogonal unitary matrices, the + 1 bases that diagonalize these classes are 
unambiguously defined and, moreover, are mutually unbiased. A maximal commut- 
ing basis of orthogonal unitary matrices is a set of + 1 sets of A^ — 1 commuting 
unitary operators (or classes) plus the identity such that these A^^ operators are 
orthogonal regarding the in-product induced by the (usual operator) trace denoted 
tr.. Let us introduce the generalised displacement operators defined as follows: 

N-l 

V^/= E7jJ'®"^^®"^'^|A:©G^)(fc|;^,j:0...Ar-l (6) 

fc=0 

It can be shown IH [7] that the A^^ operators defined in this way constitute the so- 
called generalised Pauh group. 

They are unitary with (V"/)+ = (V/)^^ = 7®°^*®'^-'^V^^*, and are traceless oper- 
ators (if we except the identity): trVl = N.6i^o.6jfl. Besides, they obey the following 
composition law: 

Making use of the composition law (IZ|), it is straightforward to check that the 

A^ classes of operators VJ*'* ^■*®''' (where / varies from to A^ — 1 and i from 1 to A^) 
contain A^ — 1 commuting operators plus the identity, and that this is also true for 
the class Vq (where / varies from to A^ — 1). Besides, the relation tr.{{V/)~^ .Vi^) = 
N.6i^i.6j^k is valid so that the A^^ V operators form a maximal commuting basis 
of unitary operators. According to the results established in the reference it is 
sufficient to find the common eigenstates of these A^ + 1 classes of V operators in 
order to determine the value of the states of the ith MUB. Before we do so, it is 
useful to introduce the t/'s operators: these operators are equal to the V's operators, 
up to a well-chosen phase. 

We also impose that f/° = Vq. 

This phase is chosen in such a way that, inside a class of commuting operators 
the t/'s operators obey an exact groupal composition law, 

Ul-Ul = Uleai. (9) 

Actually it is the requirement of the fulfillment of an exact groupal composition 
law between the f/'s operators that partially rises the ambiguity in the determination 
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of the sign of the square root factor (7^^* 1)009009)^^ prime odd dimensions our 
definition coincides with Weyl's definition It was shown in Ref.[7j that there 

exist exactly N different choices of the phases t/{/V^^*~^^®'^' that are consistent with 
the requirement that the composition law inside a subgroup of the generalised Pauli 
group is exact. We also showed that there are N different phase choices that preserve 
the composition law and lead to bases that differ only by a Galois translation in the 
label of the basis states. These choices are necessarily in one to one correspondence^. 
In a sense this is not astonishing because the concept of MUB is independent on the 
ordering of the bases. 

Moreover, as a consequence of the composition law (jH)), the following identity is 
satisfied: 

(^^(i-l)0G«10G;i-ji (^^O-l)0G;20G'2-ji ^0-1)00(^100/2) ^ (^^{i-l)0G(^leG/2)0G(/l©G^2)^| 1 

Formally we can rewrite the previous equation as follows: (7^ i)0G(a©G6)0G(aeG'')^ i _ 
(^O-i)0G(a0Ga)^ 1 ^^u-i)Qa{bQGb)^i ^^2.(O-i)0oa0Gfe)^ i ^ ^^-^^ reminisceut of the equa- 

tion although we are dealing here with half integer powers of 7^ instead of integer 
powers. 

Making use of the identity (fTUI) . we can now check by direct substitution of the 
expression that the states |e^) are common eigenstates of the ith class {i : 1...N): 

N-l -r N-l 



k'=0 V ^ q=Q 



■1 N-l 

^ AiqS)Gi)QGii-i)QGi).,eGqQGkf(ii-i)QGqeGq)\h^o \ 

7g 1G ilG ) I'^qS^Gl' 



N ,=0 



1 N-l 

Al@Gk)f((i-l)eGl@Gl))\^^_ V- ^0G(<?©GO0Gfc/^,((i-l)0G(g©GO0G(g©GO)NiuO 

7g 17g /ff 2^ 7g 17g J FgeciA-'--'-; 



^ q®Gl=0 

The states of the computational basis are eigenstates of the operators of the 0th 
class. 



4 Bell states and MUB's. 

There is a one-to-one correspondence between (generalized) Bell states and the gen- 
eralised Pauli group ( [TJ |2] see also [H] for a different approach based on additive 

^Profs. Klimov and de Guise drew my attention on the problem caused by this ambiguity 
in relation with the different phase choices for the generalised Pauli group (during the ICSSUR 
conference hold in Besancon in May 2005). 
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and multiplicative characters of the Galois field). This correspondence is a key con- 
cept for explaining several important applications of quantum information science, 
such as quantum teleportation, quantum dense coding, quantum cloning, and it 
also leads, combined to the properties of invariance of the Bell states in MUB's, to 
a solution of the Mean King's problem valid in prime power dimensions as we shall 
show in a next section. 

Following Refs.fT^. P, and |2], we can define the generalized Bell states as 
follows: 



N-l 

\B^*,n)=N-'/'Y.^G''"^^\k*)\k®Gm) (12) 

A:=0 

In this definition, we introduced the basis states \ k*) which belong to the complex 
conjugate basis of the direct basis \k). This does not make any difference when |A;) 
is the reference (computational) basis but it does when the Bell states are defined 
relatively to a basis that possesses states with complex amplitudes when they are 
expanded in the computational basis. Formally, states of the form \k*) transform 
like bra's when states \k) transform like kets. 

According to the Eq. ©, = 7jr®''"^E^=oSc?®''"^ I'^ ©G m){k\. Up to a 
global phase and a normalisation factor. Bell states and displacement {V) operators 
are thus seen to be one and the same object. Bell states possess many interesting 
properties, they are maximally entangled and form an orthonormal basis of the N"^ 
dimensional Hilbert space obtained by taking the tensor product of the dimen- 
sional Hilbert space with itself (system + ancilla). 

These Bell states possess two appealing symmetries in relation with the gener- 
alised Pauli group: 

• the Bell states are invariant under the elements of the group (the V^s operators) 
(up to a global phase): 

Vj\BU(Vjr = 7jr®°^®°"®"^'M5^*,n) (13) 

where the upper index refers to the reference basis relatively to which the Bell 
states are defined: 

N-l 

|i?^,J = iV"^/^E7r""^|er)|efe,J (14) 

1=0 

The equality shows that (up to irrelevant, global phases) the MUB's are 
eigen bases of a finite symmetry group: a set of transformations that preserves 
the Bell states (expressed in the computational basis). 
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the Bell states are permuted among each other when they are reexpressed in 
any of the MUB's: 



= 7jf "'"®""^.(7r '^^""''""¥|i?°.,ec..ec((.-i)0on)), ik-l) = 0...N - 1(15) 

The proof is straightforward and is valid in even and odd prime power dimensions 
as well (in even prime power dimensions it is made use of the identity (jlOp ): 

7V-1 
1=0 

N-1 -I N-1 1 Af-1 

^ JY-1/2 ^ y([®Gn) ^ ^,g0G/j-^,(eG(fc-l)0G'?0G'?)-j^|^O\ ^ ^ ^0Gg'0GaeGm) (^^((fc-l)0G'?'0Gg') -) f | gO^ 



1=0 g=o VN 

N-1 1 

_ Ar-1/2 ^AeGq'eGm) {iQGn) (qeGq')eGi(Ak-i)eGiq'eGq'eGqeGg)\^\^o\\^o\ 

— l^^G 7g 7g 17g J Pg/Pg'/ 

l,q,q'=0 

N-1 1 

_ Ar-1/2 AeG{q®Gi)eGm)(nQGi)&Gi{Ak-i)eG{{q®Gi)eGiq®Gi)eGeGqeGq)\^\^o\\^o \ 

- 2^ ^7g 7g 17g J Pg/PgeG*/ 

i,q,q'=q®Gi-0 

N-1 1 

^^-1/2 ^ ^iV.W.o7r^^^^^^^^'"^S-^^^^^^^^(7r^^^^^^^¥|e;)|e;e,.) 

= ^(0G«.0Gn)_^^«.-l)0G.0Gn))i|5O^^_^^^^^^^^_^^^^^^ _ ^) ^ _ ^(^g) 

This is a one-to-one mapping between the Bell states (up to global phases). 

Remark that as the Bell states are in one-to-one correspondence with the ^'s op- 
erators, the corresponding transformation rule is valid for those operators: V^(0) = 

pW.l-e-„g,^(,_i)Q^„(0, where = E£o' 7jJ'®^"^^®''"^|eL^^) (e°,| andK^(^) = 

Ylk=o 'Ig''^'^'^^'^'^^^ \^k(BGrri) {^II't '■ ^■■■N. In prime dimensions, the invariance of the 
generalised Pauli group under conjugation by any unitary matrix that maps the 
computational basis onto a MUB is a basic property of a larger group that is known 
as the Clifford group and has many applications in number theory and quantum 
computing. [111 HZ! 

Now that we derived the transformation law of the Bell states, we have at our 
disposal nearly all the tools necessary in order to derive the solution of the Mean 
King's problem in prime power dimensions. 



5 Solutions of the Mean King's problem. 
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5.1 The qubit case. 

We shall firstly give an overview of the treatment in the simplest case (qubits).jH] 
At first sight, the problem seems to be impossible to solve, because there exists 
no common eigenstate of the non-commuting observables (Tx-, Uy and (J^, and it is 
indeed impossible to discriminate between their 6 eigenstates. The solution consists 
in introducing an ancilla and of making use of the resource provided by entanglement. 
The strategy of Alice is now to add an ancilla to the spinor that the King will measure 
and to prepare initially the maximally entangled state \BQf^A,K- In virtue of the 
identity (|T6|l . this state is invariant: \Bq q) a,k=\Bqq) a,k=\Bq^ q) a,k- When the 
King measures in the kth. basis he projects the Bell state |-Bo*,o)yi,x onto \e^*)Aiice® 

\e^i)Krng,l = ^..-N -l = l,k = Q...N =2. 

Therefore Alice must, in order to save her head, be able to distinguish between 
the 6 product states \e^*)King ® l^i) Alice, k = 0...N = 2,i = 0, N — 1 = 1. According 
to our previous conventions, the indices k = 0,1 and 2 correspond to the Z, X and 
Y bases respectively. We rewrote in terms of the Bell states the solution derived 
by B-G Englert and Y. Aharonovj^ which consists of the following: after the King 
performed his measurement, Alice performs a von Neumann measurement in order 
to measure a non-degenerate observable which is diagonal in the 4-dimensional basis 
that is defined as follows: 



l*)f = 




+ \Bq^i) A,K + -Bfo)^,i^ + ^ -Sfi)^,^) 


(17) 


l^)f = 


-^i\Bofi)A,K 


+ \Bq i)a,k — \Bi q)a,k — ^ 




l^)f = 


■^i\Bofi)A,K 


~ \Bqa)a,k + -Bfo)^,^ ~ ABi,i)a,k) 




l^)f = 


-^i\Bofl)A,K 


~ \Bo^i)a,k — \Bi q)a,k + ^ 





According to our previous definitions, the four qubit Bell states are defined as 
follows: 

(18) 

i)i|o)|} 

10 
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where m,n G {0, 1}. Consequently: 



\BIo)a,k = ^{|0)I|0)| + |1)I|1)|}, \B^,,)a,k = ^{|0)I|0)| - 
\BIo)a,k = -^{|0)I|1)| + |l)f |0)f }, \BI,)a,k = -^{|0)I|1)| - I 



It is easy to check that the two last states are orthogonal to the product 
state \e^) Alice ® \^^)King^ and the two first ones to \e\) AUce ® \e\)King- For instance 
(^f koA^oii') = (^f kiA^iE") = 0- Therefore, if Ahce observes one of the two last 
(first) states and that afterwards the King reveals that he observed a state in the Z 
basis she can infer unambiguously that the Mean King observed-measured-prepared 
the state \e\) a® \e^i)K (I^IDa® \e^o)K)- Now, in virtue of the identity (jTEI), Bell states 
transform into Bell states when we pass from one of the three MUB's to another 
one: 

\Bifi)A,K = \Bqq)a,K = \Bq*q)a,K 
\Bo,l) A,K = \Bi^q)a,K = \BI,^o)a,K 

\Biq)a,k = \Bq^^)a,k = '^|-Bp,i)a,x 
\BI,)a,k = -|5fi)A,K = {-t)\B^,^,)A,K (19) 

So that the four states l"^)^ are bijectively transformed in the four states I^E')'''' 
and |\E')^(up to unobservable phase changes!). Therefore Alice can infer without 
error the values of the spins along three orthogonal directions (a rather counterin- 
tuitive result!) and consequently save her head. It is worth noting that the solution 
expressed in Eq. (fTTj) is equivalent to the one given in Refs. [SI E] excepted that now 
it is formulated in terms of Bell states. 



5.2 The Mean King's problem in prime power dimensions. 

We shall now consider the most general case (prime power dimensions) and show 
how the properties of invariance of the Bell states in MUB's also lead to a compact 
expression of the solution, valid in all prime power dimensions. |18j 

As before, Alice prepares initially the Bell state |-Bqq)a,x- The problem to solve 

can be formulated in this way: in a Hilbert space of dimension N"^, with N a prime 
power, is it possible to perform a single von Neumann measurement that allows us 
to discriminate between the product states |e")|e^) {k = 0...N — l,i = O...A^)???. If 
this problem admits a solution, then it is possible to ascertain the value of the spin 
of a spin s particle (with 2s + 1 = A^) in + 1 MUB's. 

As in the qubit case, entanglement plays a central role in this approach, as 
well as the fact that there exists a finite field with A^ elements when A^ is a prime 
power. We shall establish that a von Neumann measurement that satisfies all the 
constraints can be realised in the following basis: 

1 N-l 

\^)iiui2) - 1^ 17g )^\Bm,n)A,K) 

-'^ m,n=0 
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where 00^ represents the multiphcation of the elements of the field with A^^ el- 
ements. This field is a quadratic extension of the field with N elements and its 
elements can be represented by couples of elements of the field with N elements. 
The procedure of quadratic extension required to pass from a field with ele- 
ments to a field with A^^ elements is similar by many aspects to the procedure 
of extension of the field of real numbers that leads to the complex field. It 
is easy to check the orthonormality of this basis, in virtue of the identity 

y-N-l -,((ii,j2)00G(ni,n)) _ at2x, r. 
Z^m,n=0 I ~ "ji,0"j2,0- 

Now, if the King prepares a product of states of the computational basis \e'}*) Alice® 
\e'i)King, I = 0...N — 1 = 1, we have that 

/iTfO IpO* \ _ ±^N-l r ^^{eahQam) {eGi2QGnQGR) mo |„0* \ ^V,p^p D 
\^ {h,i2)MA^lK/ — N 2^m,n=0"rnfil 7 \-'^0,n\^lA^lK / WUCre JX 

is the remainder of (0, 1) 0g(O, 1) after division by A^. It can be shown that 
this remainder always differs from 0, as a consequence of the fact that a quadratic 
extension of a field is a field. Finally, making use of (-Bo^nle^^eJ*^) = -^^(©g"®gO 

and of the identity ©, we get that (*(i,,i2)|e['^e['j^) = -^5Q^i2eGR,i- This shows 
that, when a detector corresponding to any state that belongs to the basis |^)('j^,j2) 
fires, we can infer unambiguously the value / that was observed (prepared) by the 
King. In order to be able to infer the value of the King's observation/preparation 
in any MUB, we must have a similar relation when we reexpress the King's states 
and the |^)('ii,j2) state in the kth basis. The transformation law for the |\I') states is 
the following: 

1 N-1 

\^){ii,i2) - 1^ 17g J'-17g )^\'^{eGn®G(k-l)&Gni)*,m}){^^) 

m,n=0 

When the King prepares a product of states of the kth basis \ef*) Alice® \Gf)King, the 
preparation contains only superpositions of Bell states of the type |-Bq, „), which 
imposes that Qcn ©g {k — 1) 0g m = 0. Thus {k — 1) 0g m = n, and 

1 N-1 

hT,\fc _ (il,i2)Q&G(m,{k-l)eGm),r^k \\ 

m=0 

+contributions of weigth 
k = 1...N 

The last expression is similar to the one obtained in the computational basis. This 
proves the "invariance" of the relevant components of the states |\I'). By a com- 
putation similar to the one performed in the computational basis, we get that 
i'^{iui2)\^i*Aiice(^'iKing) = ^^*ieG(»20G(fc-i)0GK),©G«- O^ce again, this relation allows 
Alice to infer the label / of the King's observation/preparation unambiguously from 
the labels ii,i2 of the detector that fires. 
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We can understand better the invariance of the states if we note that when 
we pass from the computational basis to any of the corresponding MUB's, the Bell 
states transform in such a way that 



G 



5l \ 



In odd prime power dimensions, (7^'^"^®'^"''') 2 .(^^^'^ i)0Gm,0Gm)-j 1 



(7r^^"'¥|5L*,n')),with 



{0G"-©G(fc-l)0Gm 



[n' ) 



/(k-i) eGi\ 
V 1 / 



(22) 



This transformation is remarkable in the sense that it preserves the symplectic 
form mi©G^2 0G^i©G"^2- Indeed, as a consequence of the transformation law 
this form is the same in all MUB's: m[ ©g ^2 0g n[ Qg "^2 = "^1 ©g ^2 ©g Qg ^2- 

In odd prime dimensions, where the Galois operations reduce to the modulo N 
operations, this property was intensively studied, [ini HZl 1201 ^ ^^^^ ^ numerous 
applications. Here we see that these properties also generalise in odd prime power 
dimensions. 

It is possible to explicit the invariance of the states 1^1/) j^) ^^re elegantly, by 

expressing it directly in terms of the conserved symplectic form. This can be done 
by relabelling them as follows: 

|iT/'\0 — 1 /V^A^-l ^,(«20Gni©Gjl0Gn)/.,(m0Gn)^|i I rO \ \ 

This relabelling is bijective because 

I*')?ii,i2) = l*)°2,0Gn/Gi?) ^' remainder of (0, 1)00g(O, 1) after division 
by always differs from 0, as a consequence of the fact that a quadratic extension 
of a field is a field. 

The transformation law for the relabelled states is 



mk,^ = mi,.) (23) 

where the transformation law between primed and unprimed paris of indices obeys 
the bijective transformation law (j22|) . 

Unfortunately this elegant transformation law is not respected in general in even 
prime dimensions, as shows our example in the qubit case where one can check for 
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oneself that the four \1/ states of Ahce's basis do not exactly transform into each 
other when they are reexpressed into different MUB's. This is also the case for the 
f/'s operators (the Pauli operators in this case) which are permuted among each 
other up to multiples of minus 1 (the additive character of the field). We mentioned 
already that such phases can be compensated by a relabelling of the basis states. 
Indeed, the group structure is preserved under conjugation so that the phases of 
the U operators expressed in a different basis may be compensated by a Galois 
shift of the label of the basis states. Nevertheless this is true for subgroups of 
the full Pauli group only and it is impossible to compensate all phases. In other 
words no relabelling of the basis states is consistent with the requirement that 

^^{eGrn0Gn)^l _^^({/c-l)0Gm0G-))l_(^(rn'0Gn'))i ^-^^ ^^/^ ^^g^^^ gq. (El) for 

all couples (m, n) SIMULTANEOUSLY, by a single relabelling of the basis states. 
For instance if we reorder the spin up and down states along the X direction (which 
means that we perform a spin flip) we change the sign of two Bell states (in the 
second terms of the inequalities ()19|l ) and not only one as required. This is a simple 
example which shows that the symplectic structure is invariant in odd prime power 
dimensions only ^. 



6 Connection with the discrete Weyl and Wigner 
distributions. 



We showed that, up to a phase, generalised Bell states are in one to one corre- 
spondence with generalised displacement operators (the V ot U operators). The 
A^^ displacement opertors form a basis of the space of linear NxN operators, and 
any density matrix can be expressed as a linear combination of the displacement 
operators, which presents interesting consequences concerning tomographic applica- 
tions. As the displacement operators are diagonal in the MUB's, there is a direct 
connection between MUB's and tomography, a fact that was already recognised and 
studied in depth in the past. [312] 

Here we shall focus on the interrelations between our results and discrete quasi- 
distributions. In odd prime dimensions, the Galois and modulo A^ operations coin- 
cide, and the discrete version of Weyl's and Wigner's distributions is well-known in 
that case. [201 The identification with our results is straightforward. For instance, 
the displacement operators defined in Eqn. (28) of Ref. pli coincide with the U 
operators defined in our paper, when the dimension is odd and prime. The Weyl 
function is intimately related to the displacement operators and can thus be ex- 
pressed in terms of the U operators. The Weyl function of a linear NxN operator 
O is defined in Ref.|21J (Eqn. 82) by the relation 

^Recently, D. Gross drew our attention on the fact that this property is also true in arbitrary 
dimensions, when we work with the modulo N operations: the even dimensional case does not 
present the same guarantee regarding the symplectic invariance as the odd one. 
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W{0,a,l3) = ^-^od.N»-rnod.Nl3/rnod.N2 Y^N^l ^-rnod.NP-n.od.Nl(^Q0\O\e'^_^_^^^^ (valid ill 

odd prime dimensions). It is easy to check that Tr.Ul-O = W{0,1, {i — l).mod.Nl), 
so that the Weyl's function corresponding to the operator O is nothing else than the 
amphtude of this operator when it is expanded in the U operators basis. This is also 
true when O is a density matrix and can be generalised in a straightforward manner 
to prime power dimensions (replacing modulo operations by Galois operations). 
In this perspective, our non-standard generalisation of the displacement operators 
to prime power dimensions provides a discrete version of "a la Galois" Weyl's quasi- 
distribution. 

As we noted before there is a one to one correspondence between Bell states 
and displacement operators. This correspondence works in both directions and we 
can thus associate to the basis of the quA^^it space a basis of linear A^xA^ 

operators. As we shall now show, the amplitudes of the expansion of a linearA^xA^ 
operator O in this basis are, in odd prime dimensions, equivalent to the Wigner func- 
tion W{0,a, (3) defined in Ref. |2I]. Expressing by the symbol the operator 

obtained from 1^1/')°^ by replacing formally Alice's conjugate kets by bras, we get 

that iy(0,a,/3) = Tr.^l^ -^yO = (l/A^) ^^"^0 7"™' '^''"™' '^^'™' '^''^(C', «, /?) 
in agreement with the relation (83) in Ref. This is also true when O is a density 
matrix so that our non-standard generalisation of the ^ -^^ operators to prime 
power dimensions provides a discrete version of Wigner's quasi-distribution. 

In odd prime dimensions, it is also possible to express Wigner operators in terms 
of the parity operator Po,o and of the displacement operators, in accordance with 
the expressions (28), (61) and (66) of Ref.j^I] that can be respectively rewritten, 
according to our conventions, as follows: 

rtf„ 1 a /■ ^ \ r\ A / .Qmod.N{{i—^)&mod.Nl&mod.Nl)\\\r{i — ^)&mod.Nl 

U[a - l,fJ - [I - ij (£)mod.N I') - [imod.N " 

P{a,(3) = D{2 0w.^«,2 0„o,,;v/?)P(O,O)(V^61) and 

P(0,0) = j;jYJl^}^oD{a, (3){yQQ) where the labels V in the numbering of the 
previous expressions refers to corresponding equations in Prof. Vourdas's paper. |2T| 

The Mean King operators can be written under the form = X^^^Iq 

(7wT°®™°''^'^®'"°''^"®'"°''^^) (Tw"/^"®'"""'^'^^)^ Vl- Reexpressing the variables 
(a, /3) via the substitution a = 2 Qmod.N a ®mod.N a' and /3 = 2 Qmod.N l3 ®mod.N P' 
we get after some simplifications the relation 

iTrO Y^N~1 ( Qmod.N{'ioi®mod.N'^l3)\\ ( .Qmod.N (oi' &mod.N ) \ I: ( .Qmod.N {"^l^&mod.N P' ) \-\rP 

^(a,/3) — l^a' ,(3'=Q\lmod.N )^ K Imod.N )^ \lmod.N J^aJ 

with (5 = 2 Qmod.N OL Qmod.N «' and ^ = 2 Qmod.N P Qmod.N P' ■ 
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The righthand term of the previous equahty is in turn equivalent to the displaced 
parity operator defined through the relation (61) in Ref.[21, (V61) so that the Mean 
King operators are equal to the displaced parity operators of Refs- PU} 1^: "^^afs) ~ 
P{a,(3)i when the dimension is odd and prime. 

It is worth noting that, although the expression of the Mean King operators in 
terms of the displaced parity operator can be generalised in a straightforward manner 
to the odd prime power dimensional case, this is no longer true in even prime power 
dimensions because the parity operator is then equal to the identity, so that the P 
operators defined through the relation (\^61) coincide with the displacement or Weyl 
operators ^. Nevertheless, the Mean King operators or generalised Wigner operators 
can still be defined, also in even prime power dimensions through the relation 



1 JV-i 

= ^ E 7i^^^®^"®^^^^^"(7jr®""¥<. (24) 

m,?i=0 

In even prime power dimensions, we obtain via this relation good candidates for 
a discrete Wigner distribution that differ from the discrete Weyl distribution. For 
instance in dimension 2 we obtain that the four Wigner operators are equal to the 
4x4 Hadamard transform of the Pauli operators: 

2^('o,o) = (^0,0 + 0-1,0 + 0-0,1 + 0-1,1), 

2^E'o,i) = (^^0,0 - 0-1,0 + 0-0,1 - 0-1,1), 

2^?i,o) = (^^0,0 + 0-1,0 - f^o,i - cTi.i), 

2*('i,i) = (0-0,0 - 0-1,0 - 0-0,1 + 0-1,1), 

(with (To.o = Identity = |0)(0| + |1)(1|, ao,i = oz = |0)(0| - |1)(1|, (Ti,o = ox = 
|0)(1| + |1)(0| and (Ti,i = ay = -z|0)(l| + z|l)(0|. 

One can check that the marginals exhibit the behavior that we are in right to 
expect from well-defined Wigner distributions; for instance, ^^^(0,0) + ^("0,1) ~ 2|0)(0| 
and + ^-^ = 2|1)(1|. This qubit discrete Wigner function corresponds to the 
Wootters qubit phase-space distribution [22^ and is directly related to the tetrahedron 
qubit technique of tomography developed at NUS by Englert and coworkers. j2Sl I21| 
In prime power dimensions, the generalisation is straightforward because 

Z^/3=0 ^(a,l3) — Z^/3=0 ^a,/3=0 ^'G ) ilG )^ ^ a 

''Prof. Paz signalled this problem to me during a conversation that we had at the ICSSUR 
conference hold in Besancon in May 2005 
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^eG(a0G/3) 



) 



a 



eG(a0G/3) 
G 




We can obtain similar results in the X and Y bases by realising the two other 
possible splittings of the set of four Wigner operators into two pairs of operators. 
These splittings correspond to the concept of striation that was studied in depth in 
Refs. \22\ 123 |2n|- Each striation can be put in one to one correspondence with a 
group of iV — 1 commuting displacement operators plus the identity or, equivalently, 
with the associated MUB in which this group is diagonal. The marginal property of 
the discrete Wigner distribution can thus be rewritten in any of the + 1 MUBs, 
not only in the computational and dual bases. This is easily established making use 
of the identities ^ or 

It is worth noting that in our approach the ambiguities related to the derivation 
of a discrete phase space are implicitly evacuated from the beginning because we 
impose an a priori order to the MUBs (and thus to the associated striations). This 
particular choice is a consequence of our choice of the basis for the Galois field, of 
our identification between each element of the field and a state of the computational 
basis, of our particular phase choice for the U operators and of our numbering of 
the + 1 families of commuting displacement operators. All our choices were fully 
arbitrary, which refiects our intimate conviction that there is no prefered manner 
to order the set of MUB's or to order the states inside a given basis (up to a 
Galois translation). We systematically chose the ordering that seemed to be the 
most convenient and natural. In last resort, such a choice ought to be dictated by 
external arguments and symmetries; for instance if we identify the A^ dimensions 
with a physical parameter like a discretised position or angle, then the physical 
realisation of the states imposes a natural ordering, which was not the case in our 
approach. 

In any case, the transformation law (|23|) shows that, in odd prime power dimen- 
sions, an essential invariance is guaranteed: the phase space picture is the same (up 
to a bijective relabelling) whenever we pass from the computational basis to any of 
the A^mutually unbiased bases defined in Eq. Q. In even prime power dimensions, 
the symmetry is less strong because phases plus and minus appear that cannot be 
eliminated as we noted before, but our solution of the Mean King's problem shows 
that anyhow some regularity is still present in this case. It is worth noting another 
difference between the odd and even dimensional Wigner distributions presented 
here: the displaced parity operators (which are the odd dimensional Wigner oper- 
ators) are always factorisable in all MUB's into product of local operators. This is 
not necessarily true for the Wigner operators in even prime power (2™') dimensions. 



17 



7 Connection with Aravind's general solution. 



P.K. Aravind has shown TIT] how to generahse to prime power dimensions the solution 
of Aharonov and Englert, vahd in prime dimensions, to prime power dimensions. 
It is not our goal to compare in detail our approach and his approach but rather 
to situate the approaches relatively to each other. The expression of Aravind for 
the Mean King's basis is very general because it can be expressed in terms of the 
+ 1 mutually unbiased bases, independently on how they were derived. It is also 
essentially unique, up to a permutation of the state labels or a rearrangement in 
their relative positions. Our solution is less general because we chose a particular 
expression and ordering for the MUB's; nevertheless, it must coincide with Aravind's 
expression once this choice is made because Aravind's states can be shown to be 
orthogonal to the states \e-^* e."^) a,k with k{m) 7^ (m — l)kQ + ki, when m 7^ and 
m : l...A^ and k{fS) 7^ /cq, when m = 0. This constraint defines unambiguously the 
state that is labelled by a pair of indices k^.ki (up to a global phase). One can 
check that the states |^)° j^jj) defined in Eq. fl^ and correspond to our solution of 
the Mean King's problem fulfill exactly the same constraints, which establishes the 
convergence of both approaches. Nevertheless, it is worth noting that the additive 
character that appears in Aravind's expressions is not the additive character of the 
Galois field that we considered (pth root of unity) but the additive character of 
the modulo ring (A^ = p™th root of unity). Anyhow, unicity of the solution to 
Aravind's constraints imposes that both expressions must coincide. 



8 Conclusion 



We showed in another paperj2] how the construction of the MUB's can be de- 
rived from very primitive concepts: addition, multiplication, and duality. The whole 
structure can be derived in a self-consistent manner, making use of some well-known 
properties exhibited by finite fields. The basic intuition that finally led to the resolu- 
tion of the Mean King's problem is the recognition that the generalised Pauli group 
(sometimes called Heisenberg-Weyl group) is nothing else than a discrete Fourier 
transform performed on the elements of a Galois field. All the rest can be derived 
on the basis of this simple property. The result shows the power of entanglement 
seen as a resource: certain tasks that would be impossible if we remained confined 
to a A^- dimensional Hilbert space because they would contradict the uncertainty 
principle are possible provided we add an ancilla and exploit the counterintuitive 
properties offered by entanglement. 

Although our solution can be considered merely as a special case of Aravind's 
general solution, [TU] it exhibits interesting properties, due to our particular way 
of expressing mutually unbiased bases. It provides among others an interesting 
generalisation of the discrete Wigner function that was proposed in prime dimensions 
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in the past. [201 1^ It is an open question to know whether generahsed Wigner 
functions could be associated to arbitrary set of MUB's in the same sense that 
Aravind's solution for the Mean King's problem generalises our solution. Finally our 
approach allows us to generalise the Mean King's problem in odd non-necessarily 
prime dimensions as we show in appendix. 
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Reformulation of the Mean King's problem in odd dimensions. 



Although the modulo N operations form a field only when the dimension is 
prime, many properties that were valid when we made use of the Galois operations 
can be generalised in terms of the modulo operations. Because of this, we can 
generalise the Mean King's problem in arbitrary odd dimensions as follows: Alice 
must be able to ascertain the value of iV -|- 1 (non-necessarily) commuting (and 
possibly degenerated) observables. The observables are the displacement operators 
previously introduced but now they must be reexpressed in terms of the iVth root 
of unity (instead of p) and of the modulo operations. They can be shown by a 



20 



derivation entirely similar to the chain of equalities Eq. to be diagonal in the 
+ 1 (non necessarily mutually unbiased) bases defined as follows: 

_ ^0moci.JVg0mod.Jvfc^^((«-l)0mod.JVg0mod.iV'?)/mod.Jv2-j 1^0^ (25) 

where represents a A^th root of unity. As usually the 0th basis is the compu- 
tational basis. The solution of the generalised Mean King's problem in arbitrary 
odd dimensions is entirely similar to the case where dimensions are prime pow- 
ers, excepted that we must replace the Galois operations between A^ elements by 
the modulo A^ operations and the product (ii,i2) © &g{^,^) by 0mod.7v(^i Qmod.N 
n) ®mod.N (^2 Qmod.N '"^)- All the transformation properties of the Bell states when 
we pass from one MUB to another one are preserved and the full reasoning can be 
repeated integrally. The main difference is that the A^ + 1 eigenbases described by 
the expression (j25p are no longer mutually unbiased (we conjecture that at most 
p + 1 of them are mutually unbiased, where p is the smallest prime divider of A^). 
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